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The standard mean field slave-boson solution for the infinite-{7 Hubbard model is re- 
vised. A slightly modified version is proposed which includes properly the incoherent 
contribution of the localized states. In contrast to the standard mean field result, this 
new proposed solution defines a unique spectral function to be used in the calculation 
of local and not local quantities, and satisfies the correct thermodynamic relations. The 
same approach is applied also to the mean field approximation in terms of Hubbard 
operators. As a byproduct of this analysis, Luttinger's theorem is shown to be fulfilled 
in a natural way. 



1. Introduction 

A paradigmatic model to study systems of strongly correlated electrons is the 
Hubbard model. It contains the minimum of features necessary to describe ban- 
dlike itinerant or localized electrons depending on microscopic parameters. The 
Hamiltonian has the simple form: 

H = -Y,Y1 ^ ^ ^ ""^T ""'i ''^■^ ' ^ ^ ^ 

cr ij i 

where the first term is a one-band tight-binding Hamiltonian. Its first kinetic term 
describes the destruction of an electron on site j and the creation of an electron with 
the same spin on site i. The second term in equation (|^) is the on-site repulsion 
between electrons, which is expected to be relevant in real materials for d- and /- 
orbitals. While the limiting case U <^ t can easily be dealt with by means of per- 
turbation theory, leading to a renormalized band of quasiparticles, the intermediate 
and the strong-coupling limits J7 ~ t are much more interesting. From a qualitative 
point of view, in these regimes the system is described by a narrow coherent band 
of itinerant quasiparticles on the top of a large background of incoherent localized 
states. The metal-insulator transition is thus characterized by the disappearing of 
the coherent band. This picture is confirmed by dynamical mean field analysesEJ. 
In the strong coupling case {U ^ t) the bandwidth and the spectral weight of the 
itinerant quasiparticles scales only with the number of holes i5, so that an insulating 
state is achieved at half-filling 5 = 0. In such a situation the total weight of the 
spectral function is in its incoherent part. 
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On the analytical ground, all the possible informations about the single-particle 
properties of the system arc obtained by the knowledge of the one-electron Green's 
function G. Without losing any generality, it can be splitted in a coherent and an 
incoherent contributiona: 

G.(k,w) = — 1— — . + Gr(k,cj), (2) 

^ ' - e(k) + ^ + «0+ ^sgn[e(k) - ^ ' ^' w 

where Z is the spectral weight of itinerant states with effective dispersion e(k) 
and G'"'^ contains all the physics not described by the first term. For the above 
arguments, we expect Z (x 6. 

A meaningful quantity that can be calculated by the Green's function is the 
occupation number 

,v(k)= / ;^G,(k,c.)e— °". (3) 
J 27ri 

In a similar way as the Green's function, ^^(k) can be distinguish in a coherent part, 
with a sharp jump of height Z at the Fermi surface, plus the incoherent background 
spread over the whole Brillouin zone. 

One of the most nopuiar tools to deal with strongly correlated systems is the 
slave-boson techniquetl. Its feasibility makes it suitable to be applied to different 
models, as for instance the Kondo problem or the Anderson HamiltonianQ. In the 
finite U Hubbard model, it is applied by introducing four auxiliary bosonic fieldsEl, 
while in the t-J Hamiltonian, obtained as strong-coupling limit of the Hubbard one, 
only one boson is needed. For the purpose of this work we restrict our study to 
the simplest case of the so-called "t-model", equivalent to the infinite-C/ Hubbard 
model or to the t-J model with zero exchange coupling constant. In this case the 
strong correlation arises from the constraint of no double occupancy on each site. 

As above discussed, in the U oo limit the metal-insulator transition, at zero 
temperature, is driven by only one relevant parameter, the hole number S, related 
to the total number of electrons n hy 5 = I ~ n. In the mean-field approximation, 
which represents its simplest formulation, the slave-boson technique maps the strong 
correlated Hamiltonian JJito an effective model of non interacting fermions with 
renormalized bandwidthQ. Thus, the resulting Hamiltonian is thought to describe 
the physics of the itinerant coherent states, whereas the contribution of the localized 
states is totally disregarded at this level, and it can be restored only by higher order 
approximations. In agreement with this picture, the spectral weight associated with 
the coherent states can be shown to be equal to the doping 6 and it can be therefore 
identified with Z. If one considers the occupation number of the electrons, the mean- 
field solution of the slave-boson approach accounts only for the sharp jump at the 
Fermi energy while the incoherent background is missing (figure 1). Incoherent 
states are described in the slave-boson picture only at higher order than the mean- 
field solution. An analysis of the finite-?/ Hubbard model shows that the dynamics of 
the auxiliary bosons reconstructs the lower and upper Hubbard bands with splitting 
~ iM. However in the infinite-?/ Hubbard model, where single occupancy constraint 
is more compelling and no energy scale but the kinetic one exists, it is not clear 
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Figure 1: Sketch of the electron occupation number. The shaded area represents 
the spectral weight of the itinerant quasi-particles as described by the mean-field 
slave-boson solution. 

where the incoherent background should be located and how the total sum rule 
fulfilled. In the t-J model, where the inclusion of the exchange term identifies a 
characteristic magnetic energy, spectral weight is expected to arise on a scale J 
associated with spin excitations. However, the violation or less of the total sum 
rule is basically connected with the single occupancy constraint, independently of 
the exchange term which would just redistribute the total spectral weight. In this 
perspective, the simple i-model is as respresentative as the t-J in order to test the 
conservation of the total spectral weight. 

The aim of this paper is to revise the mean-field solution in order to take into 
account in a proper way the incoherent background. It is shown that, even in 
mean-field approximation, a rigorous implementation of the no double occupancy 
constraint leads to an analytic expression of the Green's function that describes 
on the same foot itinerant and localized states. According to its mean-field nature 
the proposed approach shares with the standard mean-field solution the common 
shortcomings due to neglecting the boson dynamics. In particular incoherent states 
in lower and upper bands can not be described. Nevertheless it provides an im- 
proved version of the mean-field slave-boson solution free of some inconsistencies of 
the standard one and preserving the correct spectral weight sum rule. The same 
procedure is applied to the "mean-field" solution of the t-model in terms of Hub- 
bard X-operators, which differs deeply from the one obtained by the slave-boson 
technique. It is shown that the inclusion of the incoherent part permits to overcome 
some intrinsic inconsistencies of the X-operator solution, as for instance about the 
validity of Luttinger's theorem. 

2. Mean field approximation in the slave-boson approach 

Let us start with a brief summary of the well-known mean field solution within 
the slave-boson formalisnm This reviewing is finalized to present all the main 
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analytic steps of such a derivation permitting to point out the critical passages that 
will be afterwards modified in order to obtain a consistent picture. 
Let us consider the Hamiltonian of the i-model: 

H = ti.jc\„Cj„ - M 51 (4) 

where c, are electronic fields operating on the reduced Hilbert space with no 
double occupied sites. 

In the slave-boson formalism the creation (destruction) operator of the real 
electron is decomposed in two operators by the usual substitution: 



c. 



- hfL (5) 
Cic — > &l/i<T, (6) 

where the h and / operators fulfil respectively bosonic and fermionic algebras. The 
constraint of no double occupancy is implemented on each site by the condition 

b\br + fl^h^+flhi = l. (7) 

By using the relations (|^)-(^), and introducing the Lagrange multiplier Ai on 
each site to enforce the constraint, the slave-boson Hamiltonian reads: 

^ = E t^Mflfo- - M E + E [^\^'^ + + /i/a - 1] ■ (8) 

In the derivation of equation (0) the relation 



b^b\flh. = flh. (9) 



has been used. 

In the mean- field approximation the local fields bi{t), b\{t), \i{t) are replaced by 
their mean-field global values: bi{t) — > {bi) = b, 6|(<) = b* , Xi{t) — > (Ai) = A. 

Note that in such an approximation both the space and time dependences of b and A 
are dropped. The resulting Hamiltonian describes an effective model of renormalized 
non interacting fermions: 

H'^^ - E b't^^fU^ - (M - A) E /I/- + aE - 1] ' (10) 

or, in Fourier space, 

H^'F - E [^'^(k) - A* + A] /L/k. + A [6^ - 1] . (11) 

k,cr 

The Green's function of the /-fermions takes the simple form 

f 1 
G^(k;a;) = ^ ^^^^^^ + ^ - A + isgn[bH{k) -^l + X]' ^^^^ 
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that represents the Green's function of purely itinerant fcrmions with bandwidth 
renormahzed by a factor and chemical potential shifted by —A. The corresponding 
zero temperature occupation number of the /-particles, 

nU^)= I |^G/(k,^)e— °-, (13) 
J 2txi 

has thus a jump from 1 to at the chemical potential. For simplicity in the following 
the zero temperature case will be always considered. Discussion and results can be 
however straightforwardly generalized at finite temperature. 
The internal energy can be expressed as function of n^(k): 

E = Y, - ^ + A] n^(k) + A [6^ - 1] . (14) 

k,cr 

The physical solution for the mean-field parameters and A is obtained by mini- 
mizing equation ([ij). Their analytical expressions are: 

A = ~^i(k)n^(k), (15) 

k,(T 

and 



62 = l-^n/(k). (16) 

k,(T 

By using these relations, equation (^^ is therefore simplified: 

£; = ^62i(k)n/(k)-/i^n/(k). (17) 

k,fT k.cr 

The total number of electrons n is easily obtained by the thermodynamic relation 
n = —dE/d^i: 

n = n/^^n^(k), (18) 

k,cr 

which identifies the number of "real" electrons with the number of the /-fcrmions, 
that is to say with the number of single occupied states. 

Equation (18) is indeed consistent with relation (^). However, some intrinsic 
inconsistencies appear just by looking at the spectral function. If one considers the 
Green's function of real electrons 

ai,3\t) = -i{TtC,„(t)c]^), (19) 

the relation between Ga and G^, within the mean field approximation, is simply: 

G,{i,r,t)^b'Gi{i,j-t). (20) 

The same relation thus holds for the spectral functions yl(k, w) — iP' {]<i^ u) , and 
for the corresponding electronic occupation number 

n(k) = &2n/(k). (21) 
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This result is in partial agreement with the physical situation depicted in figure 
1, where the contribution of the coherent states to the total number of electrons 
scales with Z = 5 = 1 — n, but it is in open contrast with Equation (]T8|). In order 
to reconcile the two descriptions one should assume two relations linking n(k) with 
n-''(k), respectively n(k) — n-^Ck) and n(k) = 6^n-''(k). The first one should be 
employed to calculate only local quantities as the total number of electrons n, as it 
has implicitly done in the expression of the internal energy in equation (17). The 
second one should be used for non local quantities, as shown by the kinetic energy 
in the same equation. It is not clear in this procedure which of them is the physical 
occupation number of the real electrons. 

The inconsistency is clearly related to the neglecting of the incoherent back- 
ground of localized states. The spectral function A(k,uj) as well as the occupation 
number n(k) of the real electrons should contain both the coherent and incoher- 
ent contributions. In the following it is shown how the previous approach can be 
modified to take into account in a simply way both the coherent and incoherent 
parts, leading to an unambiguous determination of the spectral function and of the 
related quantities n{h) and n. 



3. Modified mean-field approach 

Actually, the procedure here proposed is quite simple. The basilar considera- 
tion arises from the realizing that the above discrepancies stem from the different 
implementations of the no double occupancy constraint. Indeed, in the evaluation 
of the local quantity cl^Cia the condition has been used to identity {c\^Ci„) with 
ifiaficr)- A- similar equality however does not hold for the non local Green's function 
of the real electrons defined in equation (|l9|) . In this situation equation can not 
be applied whereas the relation (|2^) is rather fulfilled. An appropriate way to deal 
with the generic Green's function is to split it in the local and non local part: 

G^{i,j;t) = ^i{Tt£,^{t)d]^) [1 - - ^{TtcUt)cl)^^3■ (22) 

Let us now introduce the slave-boson formalism and perform the mean field 
approximation in two steps: firstly the fields bi{t), 6|(i) are replaced by their mean 
field values with respect to the time bi{t) bi, &|(i) — > &|. Note that they still 
preserve a full site dependence. As a consequence of this first approximation the 
relation (|^) can now be formulated in the dynamic version: 

5lm. = hb\fi{t)h„ = fl{t)h.. (23) 

We can now safely apply in equation (|2^ ) the spatial mean field approximation 
in the first non local part, while in the second local term the relation ( [2^ ) can 
correctly be employed. The resulting Green's function becomes: 

G,{i,j-t) = -ibHTthAt)fl) [1 - - i{TtU.{t)fl)5,j, (24) 

or, in a more compact way, 

G,{i,j-t) = -ib'{TtU{t)fl) b'] {Ttn,{t)fl)S,,. (25) 
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Figure 2: Occupation number of real electrons as calculated by equation (29). It is 



composed by a coherent part with Zn^ states and an incoherent one with (1 — Z)n 
states. 



In Fourier space, the analytic expression of the electronic Green's function takes 
the form: 

G,(k; ^) = ZGiik; c^) + (1 - Z) ^ G,(k; cu), (26) 

k 



where Z is the coherent spectral weight Z — b^. A straightforward consequence of 
equation (Kg) is: 



k 



G.(k;^) 



(27) 



so that the electronic propagator G can be fully expressed as function of G^: 

G,(k; u;) = ZGlik; u;) + {1 - Z)J2 (k; ^). (28) 

k 



From equation 
diately evaluated: 



the occupation number of the real electrons is also imme- 



rik) = Zniik) + (1 - Z) n,, 



(29) 



where ricr is the total number of electrons per spin. 

The analytic expressions of the Green's function for the /-fcrmions, and con- 
sequently of nl(k), in the present approach are just the same as in the standard 
mean field slave-boson result given by equation ( p^ ) . This can be easily checked by 
considering that the mean field Hamiltonian in term of /-operators, as defined in 
equation (|l^), is still valid in the present approach. 

In agreement with the intuitive picture, the occupation number derived in the 
present paper contains two contributions: one describing itinerant dispersive states 
with total spectral weight Z and accounting for the sharp jump at the Fermi energy; 
and an incoherent part of localized states with no k-dependence and containing 
(1 — Z)n spectral weight. The resulting picture is shown in figure 2. 
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It should be stressed that equation ( p8| ) defines in an unambiguous way the 
Green's function of real electrons G, although expressed as function of . It can 
be now directly employed to calculate any quantity involving one-particle physics, 
with no more need to deal with auxiliary fields. In particular, this unique definition 
of G, and consequently of n(k), yields a consistent derivation of the total number 
of electrons and of the internal energy, overcoming the discrepancy found in the 
standard mean field slave-boson solution. 

By working directly in c-operators, the total number of electrons n and the 
internal energy E are expressed as function of ^^(k) as: 

" = E(C^k..)=5]n.(k), (30) 

k,<T k,(T 

E^Y. W'') - (4,.£k,.) = Wk) - ".(k). (31) 

k,fT k,cr 



Inserting equation (|29|) in ( pQl ) we then obtain 

= ^n^(k). (32) 



n 

k,(T 



On the other hand, as previously discussed, the number of particles can be derived 
also by the thermodynamic relation n — —dE/dfi. The internal energy, by using 
the same expression of na-{\i) as in the calculation of n, reads: 

E = J2[b'm-^^]ntik). (33) 

k,(T 

In the same way as in the standard mean field slave-boson derivation, the phys- 
ical values of and A that appear as parameters in Ti^(k) are found: 

A = -^i(k)n^(k), (34) 

k,(T 

62 = l-5]n/(k). (35) 

k,(T 

Then, deriving E with respect to fi as required by the thermodynamic relation, one 
finally obtains 

n = 5]n^(k), (36) 

just like in equation (|32|). 

This result shows the fully consistency of the analytic expression of G derived 
in the present approach. The Green's function for real electrons is univocally deter- 
mined and its employment to calculate the total number of electrons and the internal 
energy obeys the thermodynamic relation dE/dji — —n. It is also worthy to note 
that Luttinger's theorem is naturally fulfilled by using the Green's function here 
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derived. This is particularly interesting since the validity of Luttingcr's theorem 
has not been imposed as a precondition to be satisfied. This will appear even more 
relevant if compared with the results of mean field theories performed within the 
X-operator formalism, where not only the thermodynamic relation dE/df-i = — n 
is not fulfilled, but the validity or the breakdown of Luttinger's theorem is contro- 
versial. In the next section it will be shown how the inclusion of the incoherent 
background, in the same spirit of the present analysis, allows to overcome all these 
inconsistencies also in the mean field solutions based on the Hubbard X-operators. 



Moreover the more restrictive relation 



4. Hubbard operators and mean field theories 

One of the advantages of using X-operators is that they permit to work directly 
in terms of real electrons, without introducing any auxiliary field. In the case 
of the t-model here considered, the X-operators live in the reduced Hilbert space 
constructed on each site by the states |0), | t), | i)- The Hubbard X-operators can 
then be represented as projection operators: X°''^ — \a.){P\. It is easy to check that 
they follow the algebra: 

- % (Sp^Xf ± 5o.&xf^ . (37) 
ion 

= 5p^X°'^ (38) 

is obeyed. 

By using the Hubbard operator formalism, the t-model Hamiltonian can be 
written as: 

H=Y1 - M E ^r- (39) 

i^j.(T i,cr 

In similar way, the Green's function of real electrons is defined as: 

G4t,r,t)^~i{TtX^^t)Xf). (40) 

Different mean field approximations in terms of Hubbard X-operators have been 
performed in literature by using several approaches: diagrammatic studiesQu as well 
as derivations based on the equation of motiont2l. All these analyses converge to a 
unique expression of the mean field Green's function: 

G<,(k;cj) = — — p— — — -r, (41) 

uj — Zt[k) + /i — A + 2sgn[zt(k) — /i + X\ 

where, unlike in the slave-boson technique, 

^ = l-f, (42) 

and 

X = -^Y.t{k)nUk). (43) 

k,(T 
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Coherently with the sph'it of a mean field approximation, equation describes a 
system of non interacting electrons with total spectral weight (1 — n/2) and with a 
dispersive band renormalized by the same factor. 

The occupation number is directly evaluated by (^l])- In fact, equation ( |38| ) 
implies as a particular case the relation 

X^" = XfX^", (44) 

hence 

n,(k) = (l - ^) e [a* - (1 - n/2) tik) - A] , (45) 
where Q is the Heaviside function. The corresponding total number of electrons, 

^^=(l-i)E"-W> (46) 

k,cr 

shows one the peculiarities of the mean field solutions based on X-operators, namely 
the breakdown of Luttinger's theorem. In fact, according to equation (^), the half- 
filled case n — 1 corresponds to the complete filling of the electronic band and to a 
vanishing Fermi surface. On the contrary, the half-band filling case corresponds to 
n = 2/3 leading to unphysical results. For instance, all the properties and instabili- 
ties of the system due to possible nesting of the Fermi surface, as antiferromagnctic 
or charge-density-wave ordering, lie around this particular value n — 2/3, a clear 
artifact of the approximation. 

On the other hand, the same breakdown of Luttinger's theorem appears ques- 
tionable. One could as well use the above widely discussed thermodynamic relation 
dE/dfi = — n to determine the total number of electrons. However, it is easy to 
check that the expression ( |45| ) , once plugged in the internal energy 

does not reproduce the result of equation ( ^ ) . 

We found ourselves in the same controversial situation as for the slave-boson 
solution. As before, one should postulate that the electronic propagator in equation 
( ^l] ) (and related functions like spectral function or occupation number Ucik)) has 
to be used only for non local quantities as the kinetic energy, but it does not give 
any information about the total number of particles. This inconsistency was in that 
case solved by the introduction of the incoherent background. An identical result 
will be also recovered in the Hubbard operator approach. 

Let us rewrite the "correct" occupation number by adding an incoherent contri- 
bution B: 

n,{k) = (l - ^) e [m - (1 - n/2) t(k) -X]+B. (48) 
The total number of electrons is thus given by: 

n = (^1 - X e [/i - (1 - n/2) t(k) - A] + 2B. (49) 

k,cr 
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By substituting equation ( |4q ) in ([47|), the internal energy becomes: 

^ = (l - ?) E - (1 - "/2) - A] - /in. (50) 



k,o- 



The thermodynamic relation can now be employed to determine the spectral weight 
of the incoherent background B by requiring the total number of electrons to be 
equal to equation (p9|). The derivative of the internal energy in equation (|5^) with 
respect to n gives 

1^ = -« = - n - + (/* - A) |- ^ e [/. - (1 - n/2) t(k) - A] , (51) 
where the equality 

)^^6^,* V- ■■'/^j-y^^j -.J VA- 



(l - f ) i(k)^e [m - (1 - n/2) t(k) - A] = (m - A) —6 - (1 - n/2) t{k) A] 

(52) 

has been used. Equation (pll) can be further simplified to obtain 



a^ = 9TlE0[/^-(l--/2)*(k)"A], (53) 

k,(T 

whose solution is trivially: 

n = ^e[/i-(l-n/2)i(k)-A]. (54) 

k.cr 

The constant i? is thus easily obtained by equating (^9|) with (|5 

S = Y- (55) 

Note that the validity of Luttinger's theorem is now unambiguously defined and, 
just as in the slave-boson formalism, arises in a natural way from the correct no-(k) 
which includes the proper incoherent background. 

5. Summary and discussion 

A slightly modified version of the slave-boson mean field approximation has been 
proposed in this paper to account for the localized incoherent background missing 
in the standard mean field theories. In contrast to them, the resulting Green's 
function here derived can be safely used in the calculation of any, local or non local, 
quantity as the kinetic energy, the total number of electrons or the internal energy. 
It has also been shown that Luttinger's theorem is naturally fulfilled when such a 
Green's function is employed. 

This approach permits a particular compact and suitable way, although approxi- 
mate, to deal with the redistribution of the one-particle spectral weight due to strong 
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correlation effects. Broken symmetry states, at a mean field level, can be also taken 
into account in similar way. A straightforward generalization is the superconducting 
case. It is well known that no purely local s-wave Cooper pair can be established in 
a strongly correlated system, once an attractive interaction is taken into account, 
because of the forbidding of double occupancy. However, a common shortcoming of 
the standard mean field theories is to allow for a finite loc al ,s -wave order parameter 
as a consequence of the relaxing of the local constraintEEl. This spurious result 
disappears when the modified mean field approximation here proposed is used. 
Just in the same way as in the normal state, the anomalous propagator F{ij;t) = 
— z(c|-|-(i)c^l) can be splitted in a local and non local part. Following the previous 
procedure, it easy to obtain the corresponding BCS expression for the anomalous 
Green's function F: 

pfi. \ A(k) , A(k) 

Note that both the local part has the same spectral weight Z = &^ as the non local 
one because the relation can not be employed in this case, and the non local 
part have the same spectral weight. From equation (js^) it is thus clear that the 
condensate relative to local Cooper pairs, given by the order parameter J2\i ^O^'j 
is identically zero, as physically expected. 
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